Renormalization of interactions of ultracold atoms in simulated Rashba gauge fields 
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Interactions of ultracold atoms with Rashba spin-orbit coupling, currently being studied with 
simulated (artificial) gauge fields, have nontrivial ultraviolet and infrared behavior. Examining the 
ultraviolet structure of the Bethe-Salpeter equation, we show that the linear ultraviolet divergence in 
the bare interaction can be renormalized as usual in terms of low-energy scattering lengths, and that 
for both bosons and fermions ultraviolet logarithmic divergences are absent. Calculating the leading 
order effective interaction with full dependence on the spin-orbit coupling strength and the center-of- 
mass momentum of the colliding pair, we elucidate the relation between mean-field interactions and 
physical three-dimensional scattering lengths. As a consequence of infrared logarithmic divergences 
in the two-particle propagator, the effective interaction vanishes as the center-of-mass momentum 
approaches zero. 



I. INTRODUCTION 

Artificial gauge fields in neutral ultracold atomic sys- 
tems hold the prospect, both theoretically and experi- 
mentally, of realizing and exploring a wide variety of new 
physical systems [1|. Besides simulated Abelian gauge 
fields, which give rise to analogs of conventional magnetic 
fields, artificial non-Abelian gauge fields, with potential 
applications to quantum chromodynamics, have been the 
subject of intense study |2h18|. Certain non-Abelian 
gauge fields arc equivalent to Rashba-Dresselhaus spin- 
orbit interactions 19, 20] ; the first experimental realiza- 
tion of such a spin-orbit coupled ultracold system was 
reported in Ref. [2l|. Striped phases and topologically 
non-trivial states have been predicted, within mean-field 
theory, in atomic bosonic systems with a Rashba-type 
spin-orbit coupling [y, [^ . 

A Rashba spin-orbit interaction, Kp^ • cr, with k the 
Rashba coupling and pj^ = (pa;,py,0), leads to two non- 
quadratic branches in the single particle dispersion rela- 
tion [Eq. (|ni)]. As a consequence, the effective interac- 
tions between particles can exhibit nontrivial ultraviolet 
and infrared structure. This structure is in contrast to 
usual dilute three-dimensional ultracold atomic systems 
where mean-field interactions are proportional to three- 
dimensional scattering lengths. Our purpose in this pa- 
per is to elucidate the ultraviolet and infrared structure 
of the effective interactions in ultracold atomic systems 
with Rashba spin-orbit coupling, and their relation to 
physically observable scattering lengths, which is needed 
to analyze such systems both theoretically and experi- 
mentally. 

The ultraviolet and infrared problems have been con- 
sidered earlier. Yang and Sachdev ;22j analyzed a con- 
densed matter system with one branch with a dispersion 
relation similar to the lower energy branch with a Rashba 
interaction, e(p) — \{p— Kf'\/2m, obtaining both ultravi- 
olet and infrared cutoff-dependent logarithmic contribu- 
tions to the effective interaction. Similarly, Gopalakrish- 



nan et al. iil2| have recently considered one-loop correc- 
tions to the effective interaction in atomic systems within 
the lower energy branch of a Rashba interaction and also 
obtained a logarithmic dependence on the ultraviolet and 
infrared cutoffs. The arguments in Refs. [22] and 12] are 
mainly in two dimensions, but a similar argument applied 
to three dimensions with one branch of dispersion leads 
to the logarithmic dependence on the ultraviolet cutoff 
as well. The essential argument is that with a single 
branch, the two particle Green's function in the t-matrix 
equation has the ultraviolet behavior 
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where A is the ultraviolet momentum cutoff. We show 
here that with the doubly branched dispersion taken into 
account, logarithmic ultraviolet divergences, in fact, are 
absent in three dimensions so that the linear ultraviolet 
divergence in two-particle scattering can be nicely renor- 
malized away as usual, and the effective interaction de- 
pends solely on low-energy scattering lengths. Similarly, 
logarithmic ultraviolet divergences do not appear in the 
gap equation for fermions in paired states, and linear 
divergences can again be renormalized awa y in favor of 
scattering lengths, as discussed in Refs. (iSl - llSl 23]. 

Even though after renormalization the effective inter- 
action becomes free of ultraviolet divergences, the full re- 
lation of the effective interactions in systems with Rashba 
spin-orbit coupling to physically observable scattering 
lengths has interesting structure, owing in part to an in- 
frared logarithmic divergence in the two particle propaga- 
tor entering the t-matrix when the pair's center-of-mass 
momentum, q, goes to zero. In essence. 
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As a consequence, the effective two-particle interaction 
vanishes as q — >■ 0. This infrared divergence, which 
can be obtained within a single branch model, as in 
Refs. [ij, [23, is a result of an infinity of pairs of zero 
energy states with zero center-of-mass momentum. We 
construct here the full solution of the two-particle Bethe- 
Salpeter equation for the t-matrix for fermions, and an 
approximate t-matrix for bosons valid when the free 
space three-dimensional scattering lengths are small, in- 
cluding exactly the two branches of the single-particle 
dispersion relation and the infrared structure. 



II. PARTICLES WITH RASHBA INTERACTION 

We consider both bosons and fermions in three di- 
mensions with two internal states a and b, and with a 
symmetric Rashba spin-orbit interaction described by the 
Hamiltonian, 
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Here, m is the atomic mass, V is the volume of the sys- 
tem, flp annihilates an atom in internal state a with mo- 
mentum p, and 5p annihilates an atom in state b with 
momentum p; the ax and ay are the usual Pauli matrices 
between the internal states, and / is the two-by-two iden- 
tity matrix. We take the coupling k to be positive. The 
9aa, 566, and gab are the bare s-wave couplings between 
a-a, b-b, and a-b particles; for fermions, gaa, gbb are not 
present. 

The single particle part of the Hamiltonian can be di- 
agonalized by introducing operators ap and /3p by 
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where (j)i is the angle of Pi in the x-y plane, and where 
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where the upper signs are for bosons and the lower signs 
are for fermions in Eq. (|8]). The denominators v^ and 2 
in Eq. ([7]) are chosen so that the vertices in the Feynman 
rules are just the V^'-'s. 
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where (^ is the angle of p in the x-y plane; the Hamilto- 
nian becomes 
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The single-particle spectrum has two branches. 
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with the lower branch e-(p) having degenerate single- 
particle ground states on the circle p_L — k. 

The interaction part of the Hamiltonian, "Hint , in the 



III. T-MATRIX 

To obtain the low-energy effective interaction, we con- 
struct the t-matrix describing the collision of two atoms 
in the a-branch with incoming momenta q/2 + p and 
q/2 — p and outgoing momenta q/2 -f p' and q/2 — p'. 
The momentum of each particle is on the degenerate 
ground-state circle. The single particle propagators are 
l/(w — e±(p)), and characteristically, the interactions in 
the a-/3 basis are dependent on angle. The t-matrix is 
the sum of ladder diagrams (Fig. [IJ . We denote the mo- 
menta of particles in the intermediate off-shell states by 
q/2-f k and q/2— k and label angles of the momenta (^i as 
in the figure. The Bethe-Salpeter equation for the zero- 
energy vertex function V^ with incoming and outgoing 



a-a particles is 
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where V^p is the t-matrix for incoming /3-(3 particles col- 
hding and becoming a-a particles, FJ^o is the t-matrix for 
incoming a-/3 colliding and becoming a-a particles, and 
5 and 6 label the intermediate states. The Fog and F"^ 
obey similar Bethe-Salpeter equations. We note that p 
and p' dependences of F(p,p';q) are only through the 
angles (f>i. Counting powers in a term- by-term expansion 
of the t-matrix, one finds that the kernels are linearly di- 
vergent in the ultraviolet in each order in the expansion, 
since the denominators have two powers of k while the 
integral is over three dimensions. 
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FIG. 1. The scattering t-matrix for two particles in the a- 
branch. The solid lines denote particles in a-branch, and the 
dashed lines are particles in the /3-branch. The (j)i axe the 
angles of the corresponding momenta in the x-y plane. 



We first show that there is no ultraviolet logarith- 
mic divergence in the t-matrix. For fermions the 
Bethe-Salpeter equation for r552(p,p';q) implies that 

r22(p,p';q) = -r^^(p,p';q). Also, since vj^^^^^^^,^^ = 

—VI J. J, J, , the numerators of the first two terms in 

'?'l,02;05-06' 



Eq. ([To]) are equal, and depend on k only through the 
angles of q/2 ± k. The first two terms when combined 
have a linear but no logarithmic divergence. The denom- 
inator in the third term is explicitly a function of fc^ and 
kl for large k and, thus, the integral also does not contain 
an ultraviolet logarithmic divergence. The entire integral 
in Eq. PH)) has only a linear divergence. 

For bosons, we note that for large fc, (Jjq ^ (pk and 05 ^ 
0fc -|-7r, where <j)k is the angle of k in the x-y plane; hence. 
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leading order. Also, the Bethe-Salpeter equations for 
rj^J^(k,p';q) and r2o(k, p';q) are the same to leading 
order when k is large. Thus, the same mechanism as for 
fermions leads to a cancellation of the logarithmic diver- 
gences, leaving only with the linear divergence, on which 
we now focus. 

Looking at the Bethe-Salpeter equation as a pertur- 
bation series, we see that in the intermediate processes, 
to leading order in l/k"^ the coupling gf,;, appears in the 
vertices multiplied by a factor of e^*"^*" and the coupling 
Qab by e*"^'' , while the coupling gaa is not multiplied by 
a phase factor. This implies that in integration over in- 
termediate momenta, cross terms between different cou- 
plings do not lead to linear divergences due to the phase 
integrals, whereas the phase factors in terms with the 
same couplings are canceled and form geometric series. 

For bosons, we obtain an approximate effective inter- 
action by adding linearly divergent terms, which yield 
the leading-order terms in the physical scattering lengths, 
ttij, times K. The higher-order terms for bosons include 
cross terms between different scattering lengths, which 
do not suffer from ultraviolet divergences. The linear 
divergences sum to give the effective boson interaction: 
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defines the three physical three-dimensional scattering 
lengths, Qaa, abb, and Uab, in the absence of spin-orbit 
coupling in terms of the bare couplings and the ultravio- 
let cutoff A. The functions f{q/2) and g{q/2) are 
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where q = q/ k.. When the colliding particles have zero 
energy, the center of mass momentum q can be at most 
2k, which gives the restriction < q/2 < 1. The func- 
tions f(q/2) and g{q/2), which do not diverge in the ul- 
traviolet, are plotted in Fig. [2] The effective interaction 
Eq. pT|) is free of ultraviolet divergence. 

For fermions, where aaa = Obb = 0, no cross terms 
exist and the geometric series obtained by adding linearly 




FIG. 2. The center-of-mass momentum dependence of f{x) 
(solid line) and g{x) (dashed line) entering the renormalized 
interaction [Eq. (fTTj) ]. 



divergent terms is an exact solution of the Bethe-Salpeter 
equation: 
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Both functions f{q/2) and g{q/2) diverge logarithmi- 
cally at g = 0, as noted earlier in Eq. ([2]), leading 
to the logarithmic vanishing of the effective interaction 
as g — > 0, in agreement with the result obtained by 
Gopalakrishnan et al. [l2|. This infrared divergence 
arises from the existence of infinitely many pairs of zero- 
energy single-particle states with q = 0, which add up to 
give an infrared divergence in the two-particle propaga- 
tor in the t-matrix. For q ^ 0, there is only one pair of 
zero-energy states and, thus, no infrared divergence. 

The result derived here implies that, when the spin- 
orbit coupling is small (a^jK ^ 1), the effective inter- 



action at q away from is essentially given by replacing 
the bare couplings by scattering lengths: gij — )■ 47ray /m, 
but around q ^ 0, the effective interaction essentially 
depends on the center-of-mass momentum and vanishes 
with an inverse logarithm as 5 — >■ 0. For stronger spin- 
orbit coupling (larger UijK < 1), the center-of-mass de- 
pendence of the effective interaction comes into play for 
a large range of q. For bosons with even stronger spin- 
orbit coupling {aijK^l),we need to include higher-order 
terms in the t-matrix beyond those included in Eq. (fTTj) 
to obtain the effective interaction. 

The gap equation for fermion pairing of a-fermions 
with 6-fermions in a singlet channel has a structure sim- 
ilar to the first two terms in Eq. (jlOp . which, with the 
two terms combined, contain only a linear and not a loga- 
rithmic divergence [13l - ll8l . l23j . The linear divergence can 
be renormalized away in favor of the three-dimensional 
scattering length, a procedure requiring the existence of 
two dispersion branches. Studies with only one branch 
of the dispersion relation yields, as noted, logarithmic ul- 
traviolet divergences in the effective interaction[13. |22| . 

The implications of the effective interaction derived 
here for the ground state and the finite temperature prop- 
erties of many-body Bose and Fermi systems will be dis- 
cussed in a future publication. 
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